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Abstract 

Here we review a method for constructing exact eigenvalues and eigenf unctions of a 
many-particle quantum system, which is obtained by adding some nonhermitian but PT 
invariant (i.e., combined parity and time reversal invariant) interaction to the Calogero 
model. It is shown that such extended Calogero model leads to a real spectrum obeying 
generalised exclusion statistics. It is also found that the corresponding exchange statis- 
tics parameter differs from the exclusion statistics parameter and exhibits a 'reflection 
symmetry' provided the strength of the PT invariant interaction exceeds a critical value. 
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1 Introduction 



It is well known that integrable dynamical models and spin chains with long range 
interactions exhibit fractional statistics or generalised exclusion statistics (GES) 
which is believed to play an important role in many strongly correlated systems of 
condensed matter physics. The A^-i Calogero model (related to A^-i Lie algebra) is 
the simplest example of such dynamical model, containing N particles on a line and 
with Hamiltonian given by 0, |3| 

1 " d 2 J " 2 J„ 1 

"-2§ir^ + 2£(^r (L1) 

where g is the coupling constant associated with long-range interaction. One can exactly 
solve this Calogero model and find out the complete set of energy eigenvalues as 

E m> n 2) ..,n N = ^[1 + (N- l)p)+uJ2ny (1-2) 

i=i 

Here rijS are non-negative integer valued quantum numbers with rij < Tij+i and v is a 
real positive parameter which is related to g as 

g = v 2 -v. (1.3) 



It may be noted that, apart from a constant shift for all energy levels, the spectrum (|1 . 2| ) 
coincides with that of N number of free bosonic oscillators. Furthermore, one can easily 
remove the above mentioned constant shift for all energy levels and express (|1.2j ) exactly 
in the form of energy eigenvalues for free oscillators: i? nij „ 2j ... jnjv = 4r+w 12j=i where 
fij = rij + v(j — 1) are quasi-excitation numbers. However it is evident that these fijS 
are no longer integers and they satisfy a modified selection rule given by fij+i — n,j > v, 
which restricts the difference between the quasi-excitation numbers to be at least v 
apart. As a consequence, the Calogero model ( |1.1| ) provides a microscopic realisation 
for fractional statistics with v representing the corresponding GES parameter ||, |[ || 0. 

Recently, theoretical investigations on different nonhermitian Hamiltonians have re- 
ceived a major boost because many such systems, whenever they are invariant un- 
der combined parity and time reversal (PT) symmetry, lead to real energy eigenvalues 
[§, |, [T(], |LT| . This seems to suggest that the condition of hermiticity on a Hamiltonian 
can be replaced by the weaker condition of PT symmetry to ensure that the corre- 
sponding eigenvalues would be real ones. However, till now this is merely a conjecture 
supported by several examples. Moreover, in almost all of these examples, the Hamil- 
tonians of only one particle in one space dimension have been considered. Therefore, 
it should be interesting to test this conjecture for the cases of nonhermitian iV-particle 
Hamiltonians in one dimension which remain invariant under the PT transformation 

[0: 

% -> -i, Xj -> -Xj, Pj -> Pj , (1.4) 
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where j 6 [1, 2, ■ • • , N], and Xj (pj = —i-£r) denotes the coordinate (momentum) opera- 
tor of the j-th particle. In particular, one may construct an extension of Calogero model 
by adding to it some nonhermitian but PT invariant interaction, and enquire whether 
such extended model would lead to real spectrum. 

The purpose of the present article is to review the progress 0, [12] on the above 
mentioned problem for some special cases, where the PT invariant extension of Calogero 
model can be solved exactly. In Sec. 2 of this article we consider such a PT invariant 
extension of Aw-i Calogero model and show that, within a certain range of the related 
parameters, this extended Calogero model yields real eigenvalues. Next, in Sec. 3, we 
explore the connection of these real eigenvalues with fractional statistics. Section 4 is 
the concluding section. 

2 Exact solution of an extended Calogero model 

Let us consider a Hamiltonian of the form @ 

n = H + 5Y,^—^-, (2-1) 

jjLk X 3 X k ° x j 

where H is given by eqn. fll.lj) and 5 is a real parameter. It may be observed that though 
the Hamiltonian (|2.1| ) violates hermiticity property due to the presence of momentum 
dependent term like 5J2j^k X - Xk ^ remam s invariant under the combined PT trans- 
formation ( |1.4j ). Next we recall that, Aw-i and Calogero models as well as their 
distinguishable variants have been solved recently by mapping them to a system of free 



oscillators [13|, [TJ|, [T^, |16j . With the aim of solving the PT invariant extension ( |2.1| ) of 
t4at_i Calogero model through a similar produce, we assume that (justification for this 
assumption will be given later) the ground state wave function for this extended model 
is given by 

= e~ U=i x i J] {xj - x k y, (2.2) 
j<k 

where v is a real positive number which is related to the coupling constants g and 6 as 

g = v 2 - u (l + 25) . (2.3) 

Now if we use the expression (|2.2|) for a similarity transformation to the Hamiltonian 
O), it reduces to an 'effective Hamiltonian' of the form 



W = tp- l Hifj gr = S~ + uS 3 + E gr , (2.4) 
where the Lassalle operator (S~) and Euler operator (S 3 ) are given by 
1 N B 2 18 N 3 

s~ = ~Ep -iy-s) Y.^—^ri s 3 = ^^-, (2.5) 

/ J= l OXj j_^ k Xj Xk OXj - =l OXj 



.3 



and 

E gr =^[l + (N-l)(u-5)\. (2.6) 



It is easy to see that the Lassalle operator and Euler operator, as defined in eqn.( [2.5|) , sat- 
isfy the simple commutation relation: [S 3 , S~] = —2S~. Using therefore the well known 
Baker-Hausdorff transformation we can remove the S~ part of the effective Hamiltonian 
Ti! and through some additional similarity transformations reduce it finally to the free 
oscillator model [[/]] 

i N pa ,,2 JV at- 

H free = 5- 1 (H> - E gr ) 5 = -E^ + yEfy. (2-7) 

Z 3=1 " J J ~ 3=1 ~ 



where S = e ^ s ~e^ v \^^ and V 2 = E N - - 



N 

i*3 and V 2 = £ J=1 ^ 
Due to similarity transformations in ( |2.4p and ( |2.7| ), one may naively think that 
the eigenfunctions of the extended Calogero model (|2.1| ) can be obtained from those 
of free oscillators as: xfj ni , n2l ... tnN = ip gr S IlljLi e ~* x iH nj (xj)\, where njs are arbitrary 
non- negative integers and H nj (x 3 ) denotes the Hermite polynomial of order rij. However 
it is easy to check that, similar to the case of A?v-i Calogero model [[y|, the action of 
S on free oscillator eigenfunctions leads to a singularity unless they are symmetrised 
with respect to all coordinates. Therefore, nonsingular eigenfunctions of the extended 



Calogero model ( |2.1| ) can be obtained from the eigenfunctions of free oscillators as 

i) ni ,m,:,n N = *Pgr SA+ I JJ e"^fT n (xj) \ , (2.8) 



where A + completely symmetrises all coordinates and thus projects the distinguishable 
many-particle wave functions to the bosonic part of the Hilbert space. Evidently, the 
eigenfunctions ( |2.8| ) will be mutually independent if the excitation numbers njs obey 
the bosonic selection rule: > rij. Thus, in spite of the fact that the interacting 

Hamiltonian ( |2.1| ) is convertible to the free oscillator model, the need for symmetriza- 
tion shows that the many-particle correlation is in fact inherent in this model. The 
eigenvalues of the Hamiltonian ( |2.1| ) corresponding to the states ( |2.8j ) will naturally be 
given by 

N Nu N 

E ni>n2j ... >nN = E gr + uY,n j = — [l + {N-l){v-5)\+uY J n r (2.9) 

3=1 Z 3=1 



Since § and v are real parameters, the energy eigenvalues ( [2.9|) are also real ones. Thus 
we interesting find that the nonhermitian PT invariant Hamiltonian ( |2.1p yields a real 
spectrum. Furthermore, it is evident that for all rij = 0, the energy £ , ni n2j ... njv attains 
its minimum value E gr . At the same time, as can be easily worked out from eqn. 



the corresponding eigenfunction reduces to ijj gr ( |2.2| ). This proves that ip gr is indeed the 
ground state wave function for Hamiltonian (|2.1|) with eigenvalue E gr . 
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It may be observed that the eigenfunctions ( |2.8| ) pick up a phase factor (— l) v under 
the exchange of any two particles. Therefore, v represents the exchange statistics pa- 
rameter for the extended Calogero model ( |2.1| ). By solving the quadratic eqn.( |2.3| ), one 



can explicitly write down v as a function of g and 8 as 



v={5+\)±^g+{5+ l -y. (2.10) 

For the purpose of obtaining real eigenvalues ( |2.9|) as well as nonsingular eigenfunctions 
( |2.8| ) at the limit Xi — > Xj, we have assumed at the beginning of this section that v 
is a real positive parameter. This assumption leads to a restriction on the allowed 
values of the coupling constants g and 8 in the following way. First of all, for the 
case g < — (8 + |) 2 , eqn.( [2.10"| ) yields two imaginary solutions. Secondly, for the case 
8 < — |, > g > — (5+|) 2 eqn. (|2.101 ) yields two real but negative solutions. Inequalities 



corresponding to these two cases represent two forbidden regions of (8, g) plane which 
are excluded from our analysis. 

For the case g > with arbitrary value of 8, one gets a real positive and a real 
negative solution from eqn.( p.lO| ). The real positive solution evidently leads to physically 
acceptable set of eigenfunctions and corresponding eigenvalues within this allowed region 
of (8,g) plane. Finally we consider the parameter range 8 > — ~, > g > — (8 + |) 2 , 
for which eqn. (|2.10|) yields two real positive solutions. It is easy to see that these two 
real positive solutions are related to each other through a 'reflection symmetry' given 
by v — > 1 + 28 — v. Consequently, for each point on the (8, g) plane within this allowed 
parameter range, one obtains two different values of the exchange statistics parameter 
leading to two distinct sets of physically acceptable eigenfunctions and eigenvalues. 
Thus we curiously find that a kind of 'phase transition' occurs at the line 8 = — | on 
the (8, g) plane. For the case 8 > — ~, exchange statistics parameter shows the reflection 
symmetry when g is chosen within an interval — (| + 5) 2 < g < 0. On the other hand 
for the case 8 < — |, such reflection symmetry is lost for any possible value of g. 

We have seen in this section that, similar to the case of An_\ Calogero model, the 



extended model ( |2.1| ) can also be solved by mapping it to a system of free harmonic 
oscillators. So it is natural to enquire whether this extended model is directly related to 
the An -i Calogero model through some similarity transformation. Investigating along 
this line ||12j|, we find that 



where Y = Y[j<k( x j~ x k) S , and H' denotes the Hamiltonian of A^_x Calogero model with 
'renormalised' coupling constant given by g' = g + 8(l + 5). Due to the existence of such 
similarity transformation, one may expect that the Hamiltonians 7i and H' always lead 
to exactly same eigenvalues. However it should be noted that, within a parameter range 
given by 8 > 0, g > —8(1 + 8), there exists a positive solution of eqn. (|2.3| ) satisfying the 
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condition v — 5 < 0. Therefore, we can not get any lower bound for the corresponding 
energy eigenvalues (|2.9|) at iV — > oo limit. On the other hand, the energy eigenvalues 
( |1 . 2| ) of Ayv_i Calogero model are certainly bounded from below for all possible choice 
of N and g. So there exists a parameter range within which the spectrum of extended 
Calogero model differs qualitatively from the spectrum of the original Calogero model. 
To explain this rather unexpected result, we first observe that the renormalised coupling 
constant g' would be a positive quantity within the above mentioned parameter range. 
Consequently, the corresponding exclusion statistics parameter v', which is obtained by 
solving eqn. (|1.3| ), has one positive and one negative solution. One usually throws away 
this negative solution of z/, since the corresponding eigenfunctions become singular at 
the limit Xj — > Xk- However, by using the relation ( 2.111) , such singular eigenfunctions 
(denoted by i/j'(xi, x 2 , • • ■ , m ay now be used to construct the eigenfunctions of 
extended Calogero model (denoted by ip(xi, x 2 , ■ ■ ■ xn)) as 

ip(x 1 ,x 2r ■ ■ ,x N ) = Yi(xj ~ x k ) s ip'(x 1 ,x 2 , ■ ■ ■ ,x N ) . (2.12) 

j<k 

It can be easily checked that, due to the existence of the factor Ylj<k( x j ~ x k) S , the r.h.s. 
of the above equation becomes nonsingular at the limit Xj — > Xk • Thus we curiously find 
that singular eigenfunctions of H' can be used to generate nonsingular eigenfunctions 
of 7i through the relation ( 2.12|) . This shows that the similarity transformation (EOT ) 
is a subtle one and, within a certain parameter range, eigenvalues of extended Calogero 
model will match with those of Calogero model (having renormalised coupling constant) 
only if the corresponding unphysical eigenfunctions are taken into account. 



3 Connection with fractional statistics 

We have mentioned in Sec.l that GES can be realised microscopically in An_\ Calogero 
model with hermitian Hamiltonian. The GES parameter for this Calogero model is a 
measure of 'level repulsion' of the quantum numbers generalising the Pauli exclusion 
principle. Now for exploring GES in the case of PT invariant model ( |2.1|) , we observe 
that eqn. (|2.9|) can be rewritten exactly in the form of energy spectrum for iV free 
oscillators as 

E ni , m -n N = -z- + u E % ' 
A 3=1 

where 



n 



j 



nj + (v-S)(j-l) (3.2) 



are quasi-excitation numbers. However, from eqn.( |3.2| ) it is evident that such quasi- 
excitation numbers are no longer integers and satisfy a modified selection rule: fij+i — 
fij > v — 5 . Since the minimum difference between two consecutive UjS is given by 

u = v-S (3.3) 
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the spectrum of extended Aw-i Calogero model ( |2.1|) satisfies GES with parameter v 
[0]. Several comments about this GES parameter are in order. It may be to noted that 
for 5^0, the GES parameter v is different from the power index v, which is responsible 
for the symmetry of the wave function. Therefore we may interestingly conclude that 
unlike Calogero model, the exclusion statistics for model ( |2.1| ) differs from its exchange 
statistics. Furthermore it is already noticed that, on a region of (5, g) plane satisfying 
the inequalities 5 > 0, g > —5(1 + 5), there exists a positive solution of eqn.( p.3|) 
which yields a negative value of v. For this case, however, one does not get well defined 
thermodynamic relations at iV — > oo limit and, therefore, can not interpret v as the 
GES parameter. 

By using eqn. fl2.3|) and (|3.3| ), we find the relation 

v 2 -i> = g + 5(5 + 1), (3.4) 

which clearly describes a parabolic curve in the coupling constant plane (5, g) for any 
fixed value of v. As a consequence of this, the competing effect of the independent 
coupling constants g and 5 can make the GES feature of ( |2.1| ) much richer in comparison 
with the Calogero model. For example, while bosonic (fermionic) excitations in the 
Calogero model occur only in the absence of long-range interaction, the quasi-excitations 
in can behave as pure bosons (fermions) even in the presence of both the long-range 
interactions satisfying the constraint v(5, g) = (v(5, g) = 1). Both of these constraints 
lead to the same parabolic curve g = —5(1 + 5). A family of such parabolas with shifted 
apex points are generated for other values of v and the lowest apex point is attained at 
v — J, where the quasi-excitations would behave as semions. 

4 Conclusion 

Here we construct a many-particle quantum system (|2.1|) by adding some nonhermitian 
but combined parity and time reversal (PT) invariant interaction to the A^^\ Calogero 
model. By using appropriate similarity transformations, we are able to map this ex- 
tended Calogero model to a set of free harmonic oscillators and solve this model exactly. 
It turns out that this many-particle system with nonhermitian Hamiltonian yields a 
real spectrum. This fact supports the conjecture that the condition of hermiticity on a 
Hamiltonian can be replaced by the weaker condition of PT symmetry to ensure that 
the corresponding eigenvalues would be real ones. It is also found that the spectrum 
of extended Calogero model obeys a selection rule which leads to generalised exclusion 
statistics (GES). 

However, this extended Calogero model exhibits some remarkable properties which 
are absent in the case of usual Calogero model. For example, we curiously find that 
the GES parameter for this extended Calogero model differs from the corresponding ex- 
change statistics parameter. Moreover a 'reflection symmetry' of the exchange statistics 
parameter, which is known to exist for An_\ Calogero model, can be found in the case 
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of extended model only if the strength of PT invariant interaction exceeds a critical 
value. 

Finally we note that, it is possible to obtain another exactly solvable many-particle 
quantum system by adding some nonhermitian but PT invariant interactions to the 
Bn Calogero model (associated with Bn Lie algebra) Jl2|]. Such a PT invariant model 
also leads to real spectrum with properties quite similar to the case of extended A^-i 
Calogero model. 
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